Spin-Droplet State of an Interacting 2D Electron System 
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We report thermodynamic magnetization measurements of two-dimensional electrons in several 
high mobility Si metal-oxide-semiconductor field-effect transistors. We provide evidence for an 
easily polarizable electron state in a wide density range from insulating to deep into the metallic 
phase. The temperature and magnetic field dependence of the magnetization is consistent with the 
formation of large-spin droplets in the insulating phase. These droplets melt in the metallic phase 
with increasing density and temperature, although they survive up to large densities. 
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Magnetic ordering of a low-density electron system 
is determined by the interplay between the electronic 
Coulomb interaction and Pauli principle. As the density 
decreases, the ratio between the interaction and Fermi 
energies increases, pushing the system towards a ferro- 
magnetic instability. 

In the Hartree-Fock approximation, the Bloch instabil- 
ity, a first-order transition from unpolarized to fully po- 
larized state, happens at an unrealistically small r s ~ 2. 
In the opposite limit of short-range interaction the Stoncr 
instability, a second-order phase transition characterized 
by divergent spin susceptibility, is expected. The hierar- 
chy of these transitions is discussed in Ref. [l| within the 
RPA approximation. Numerical simulations for a clean 
single-valley two-dimensional electron system (2DES) Q 
predict a Bloch instability at r s w 25 followed by Wigner 
crystallization Q at r s « 37. 

However, at very low densities a realistic system cannot 
be treated as a clean one: even small potential fluctua- 
tions due to inevitably present disorder become domi- 
nant and lead to Anderson localization. At higher densi- 
ties intricate interplay between disorder and interactions 
manifests itself as a metal-insulator transition (MIT) at 
some density n c Q ■ Experimental observations [H and 
theoretical arguments @, Q suggest that a 2DES be- 
comes strongly nonuniform at densities lower than n c ; for 
n < n c a 2DES can be considered as consisting of weakly 
coupled disordered quantum dots. Disorder is also ex- 
pected to drive a 2D system further towards ferromag- 
netic instability 0, Q . In particular, a disordered quan- 
tum dot is predicted to have a finite spin in the ground 
state, a phenomenon analogous to the Stoner instabil- 
ity 13, 1 1 1 - Experiments on quantum dots 12 1 in GaAs 
indeed found spontaneous spin polarization at r s ~ 7.6, 
much smaller than the expected value for a clean sys- 
tem 0. 

Coulomb interactions lead to renormalization of the 
Fermi-liquid constants, notably the density of states and 
the effective g-factor, g = go/(l + Fq) w ith Stoner 
instability expected at Fq = — 1. Negative Fq can 



indeed be deduced for different 2DESs from measure- 
ments of Shubnikov-de Haas (ShdH) oscillations [13l.ll4| 



and the temperature-dependent resistivity 

See- 
ing analysis of the magnetoresistance led the authors of 

Refs. 0,[l{| to suggest a quantum phase transition into a 

ferroma gne tic state at n n . This conjecture was contested 



in Refs. [2fJ and [13|, 3, 2(Jl2l| on the basis of thermody- 
namic and ShdH measurements, respectively. It should 
be emphasized that all measurements have been done at 
relatively high magnetic fields, at which Zeeman splitting 
exceeds the temperature (q* jjlrB > fcgT); apparently for 
this reason no anticipated [22h24| divergency of the 2D 
spin susceptibility with decreasing temperature has been 
revealed experimentally until now. 

We report an observation of a spin-droplet state in low 
disorder Si metal-oxide-semiconductor field-effect tran- 
sistors (MOSFETs) on the basis of thermodynamic mag- 
netization measurements. Wc used a recharging tech- 
nique suggested in Ref. [H| for energy spectrum recon- 
struction and developed for magnetization measurements 
in Ref. (20[. In this technique, the recharging current 
between the gate and the 2DES flows in response to a 
change of the 2DES chemical potential fi caused by mod- 
ulation of the in-plane magnetic field at a constant gate 
voltage [H] 
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where n is the electron density, cq = e/Aird is the geo- 
metrical capacitance per unit area, and the capacitance 
c = (c^ 1 +e 2 dn/dfi)~ 1 includes compressibility contribu- 
tion. Note that, as long as the geometrical capacitance is 
magnetic field independent, the recharging current is pro- 
portional to dfi/dB, and therefore the recharging method 
is distinct from the magneto-capacitance one. We used Si 
MOSFETs with w 200nm gate oxide, much thicker than 
the 2DES « 4nm and thus Co is set by the gate oxide and 
is almost magnetic field independent 26]. 

Importantly, Eq. ([T]) holds even when the capacitance 
c acquires an imaginary part due to contact and channel 



resistances [26[. This allowed us to extend the measure- 
ments deep into the insulating phase, with the only con- 
straint being an ability to accurately measure recharg- 
ing current I and capacitance c at the frequency u) of 
magnetic field modulation. In practice, the sample resis- 
tance was required to be below ss 100 Mfl, which hap- 
pened, e.g., at less than half of the critical density n c 
at 4 K. The recharging current per unit area is then 
given by J = (—iu>c/e){dn/dB)8B, where SB is the 
modulation depth. By virtue of the Maxwell relation 
du/dB = —dM/dn, d\x/dB can be expressed as the 
derivative of the magnetization M per unit area with 
respect to the density. 

We measured dM/dn in several high-mobility (100 VSi 
MOSFET structures similar to those used in Refs. fl3ll20l . 
27| . Such a 2DES possesses a twofold valley degeneracy 
in addition to its spin degeneracy. The in-plane magnetic 
field was modulated at frequency w/2-7r = 4— 12 Hz with 
amplitudes up to 40 niT. Measurements were performed 
over a wide range of temperatures (0.4 — 20K), and in 
magnetic fields up to 9 T. 

To get an insight into the expected behavior of dM jdn 
as a function of external parameters, let us briefly review 
two limiting cases (i) of very high and (ii) of very low 
densities; in both cases electrons can be viewed as non- 
interacting. 

(i) For high densities, deep in the metallic regime, the 
ratio between interaction and the Fermi energies is small. 
One expects to get the density-independent Pauli spin 
susceptibility at a magnetic field below the field of full 
spin polarization 2n/ '(yguB)-, where v is the density of 
states, and as a result dM/dn = 0. 

(ii) For very low densities, each electron is localized in 
its own potential well. Interactions between electrons are 
small; hence, the electron spins should be polarized inde- 
pendently: dM/dn = ub tanh(b) , where b = nsB/ksT 
is a normalized magnetic field [281 ] . 

In Fig.[lja) we present the magnetic field dependence 
of dM/dn measured deep in the insulating phase at 
n = 5 X 10 10 cm~ 2 (n c « 8.5 x 10 10 cm -2 for this sample). 
The results are shown after subtraction of the diamag- 
netic contribution [2!| ; the subtraction does not affect 
the low-temperature results in any significant way. For 
temperatures above « 4.2 K, dM/dn is consistent with 
the expected dependence dM/dn = /istanh(6) for indi- 
vidual spins. However, as temperature decreases, (i) the 
low- field slope of dM/dn vs. b becomes much steeper 
than the one expected for an independent spin 1/2; (ii) 
dM/dn vs b becomes nonlinear; and (iii) at low tem- 
peratures dM/dn reaches a maximum at b w b* [30l |; 
this maximum significantly exceeds [Ib- The fact that 
dM/dn > ub means that an electron added to the sys- 
tem not only aligns its spin with the field but also pro- 
motes spin alignment of neighboring electrons. This is 
the "smoking gun" evidence for the ferromagnetic inter- 
action between spins. Indeed, all these observations can 
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FIG. 1. Panel (a): dM/dn vs normalized magnetic field 
b = (iBB/kBT at n = 0.5 x 10 11 cm -2 . We subtracted the dia- 
magnetic contribution estimated from the high-temperature 
data, ~ 0.04/is per tesla [2^|. Panel (b): the same as panel 
(a) at n — 1.5 x 10 11 cm~ 2 ; the subtracted diamagnetic con- 
tribution is m 0.035/is per tesla. Dashed lines in (a) and (b) 
show dM/dn for localized spins 1/2. Panel (c) dM/dn vs b at 
different densities at T = 1.8K. Note that dM/dn(b) becomes 
nonlinear at the density and temperature independent 6*. 

be simulated even within the mean field approximation, 
if one assumes that ferromagnetic interaction constant 
grows with density, see the Supplemental Material. 

At higher densities, well inside the metallic phase, e.g. 
at n — 1.5 x 10 11 cm -2 shown in Fig.QJb), low temper- 
ature dM/dn changes sign. Note that negative dM/dn 
is expected in the metallic phase, since increase in den- 
sity reduces interaction and therefore polarization of the 
2DES. Thermal fluctuations also suppress magnetic or- 
dering; therefore dM /dn becomes less negative with tem- 
perature and, at temperatures exceeding the Fermi en- 
ergy (about 10K at n = 1.5 x 10 11 cm -2 ) approaches the 
dependence expected for noninteracting electrons. 

Most importantly, the characteristic normalized mag- 
netic field b* is almost constant over a wide range of 
densities [See Fig.QJc)], from ~ n c /2 to ~ 3n c , thus 
excluding any possibility of a density-driven quantum 
phase transition into a homogenous ferromagnetic state. 
Rather it is reminiscent of the behavior of a large spin 
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FIG. 2. Phase diagram of dx/dn represented in colors for 
each temperature and density, in units of ([ib/T). Dashed 
blue line represents the density, no, at which dx/dn is zero. 
Black ♦ - depicts the MIT critical density, n c , known the from 
transport measurement. 



J = (1/2) (1/6*) - 2 system. 

Our results cannot be attributed to localized spins 
solely, whose interaction is known to be antiferromag- 
netic 



31] . A minimal model that would explain them is 
a two-phase state, consisting of electron droplets with a 
typical spin of 2 and itinerant electrons. As the density 
increases in the insulating phase, electrons join already 
existing droplets or populate new ones, which leads to a 
positive dM/dn. For even higher densities these droplets 
start to coexist with itinerant electrons; the addition of 
an electron to the system increases screening and there- 
fore depopulates the droplets, which leads to negative 
dM/dn of the same order of magnitude and the same 
characteristic magnetic field scale b* as in the insulating 
phase. 

The density at which dM/dn changes sign lies in the 
metallic phase and is temperature dependent; as T — > 0, 
it approaches the MIT critical density, n c . In order to 
show this, let us focus on the low-field slope of dM/dn 
in Fig. [2 dx/dn = d 2 M / dndB. In Fig. [21 we plot a color 
map of dx/dn for different temperatures and densities. 
The density no at which x is maximal (dx/dn = 0) ex- 
trapolates linearly to n c with decreasing temperature. 
The coincidence of the magnetic and transport critical 
density values, no and n c , at T = suggests that the two 
phenomena: the formation of easily magnetized droplets 
and the MIT are intimately interrelated. 

One would expect Curie (1/T) temperature depen- 
dence for the spin susceptibility of a droplet. In contrast, 
we found the temperature dependence of \dx/dn\ to be 
closer to 1/T 2 both in the insulating phase and in the 
metallic phase away from the dx/dn = region; see in- 
set in Fig. [31 This indicates that the number of droplets 
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FIG. 3. dx/dn for different temperatures. Solid line rep- 
resents dx/dn extracted from the Sh-dH results [l3T| . Inset 
shows temperature dependence of \dx/dn\ for two densities: 
(M)n = 5x 10 10 cm~ 2 , and (»)n = 1.5 x 10 n cm~ 2 with data 
points in the vicinity of dx/dn = omitted; the solid line is 
\d X /dn\ oc T~ 2 

increases as cx 1 /T in the temperature range of our mea- 
surements. 

In Fig. [3l we present cross sections of the data of Fig. [2] 
at several temperatures. Deep in the insulating phase 
(low densities), dx/dn is positive and for low enough 
temperatures far exceeds the Curie value for indepen- 
dent spins. As density increases, dx/dn drops, changes 
sign, reaches a minimum, and eventually goes to zero, as 
expected for Pauli susceptibility of noninteracting elec- 
trons, which is density independent. As temperature de- 
creases, the transition from low-density positive dx/dn 
to the high-density negative one becomes steeper, indi- 
cating development of a sharp cusp in x m the vicinity 
of n c at zero temperature. This behavior resembles the 
sharp drop of dR/dT with density in the vicinity of the 
MIT. 

Note that transport measurements, both Shubnikov- 
de Haas [HI (Sh-dH) and in-plane magnetoresistance 18, 
Hil , also found decrease in x with density in the metallic 
phase. However, the magnitude of dx/dn extracted from 
these measurements is much smaller than that measured 
with the recharging method, sec solid line in Fig. [3] This 
discrepancy far exceeds any possible uncertainty in our 
experiment. We believe it is a result of the fundamen- 
tal difference between the physical quantities provided 
by thermodynamic and transport methods: the thermo- 
dynamic method yields an average magnetization of all 
the electrons that are capable of recharging within the 
80 ms field modulation period. In contrast, the transport 
is influenced mostly by delocalized electrons having the 
picosecond-scale relaxation time. 

In conclusion, we presented experimental evidence for 
the existence of spin droplets in high-mobility 2D elec- 
tron layers in Si MOSFET samples on both sides of the 
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MIT. The absence of similar behavior in low- mobility 
samples, for which the MIT is much less pronounced and 
occurs at much higher critical density n c « 3 x 10 11 cm" 2 , 
emphasizes the importance of interactions for the spin- 
droplet formation. A typical total spin J = 1/26* « 2 of 
a droplet, estimated above coincides with the value pre- 
dicted within the theory [32j for quite a realistic interac- 
tion constant value Fq = —0.5, typical of this regime fl7j ]. 
The minimal size of such a droplet can be estimated as 
yj2J/n. This size, « 100 nm for n = 5 x 10 10 cm~ 2 , 
is comparable with the gate oxide thickness of 200 nm, 
which screens the potential fluctuations and, hence, sets 
their spatial scale. Our results are in qualitative agree- 
ment with numerical calculations which predict enhance- 
ment of spin polarizability by disorder [8[ . The observa- 
tions of spin droplets in the metallic phase are in line 
with the expectation of stochastically formed multispin 
fluctuations in Fermi-liquid [32| , in the insulating phase 
with the expectations for spontaneous magnetization of 
a disordered quantum dot 13, 11 \ ■ 
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FIG. 4. dM/dn calculated within the mean field approxima- 
tion. We assumed T c oc n k ; we show two sets of curves, solid 
lines for k — 2, and dotted lines for k = 1. Note similarity 
with the data presented in Fig. la of the main text. 
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Our results suggest a picture of a nonuniform state 
in which easily polarizable electron droplets coexist with 
Fermi liquid over a wide density range. The spatial or- 
der (if any) of electrons within spin droplets remains un- 
known; these might be e.g. Wigner crystallites, or even 
more sophisticatedly ordered droplet phases [33j . How- 
ever, one needs to explain why these two coexisting sub- 
systems interact so weakly. Indeed, b* does not manifest 
itself in transport measurements, and there is no satura- 
tion of phase breaking time at low temperatures due to 
scattering by the spin droplets 17 1 , as expected Q . A 
plausible explanation for this would be the large size and 
collective nature of the droplets, which suppress electron- 
droplet scattering with spin flip. 



The existence of these droplets can be checked in sev- 
eral ways; the most direct one would be magnetic force 
microscopy, somewhat analo gou s to measurements of lo- 
cal compressibility in 2DES [35j . Polarization of photo- 
luminescence signal serves as a local probe for the spin 
polarization. The temperature dependence of ESR signal 
may provide information on the typical magnetic moment 
and the signal width on the spin relaxation rates. Knight 
shift measured on 29 Si can be used for the same purpose. 
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The simplest framework to include ferromagnetic inter- 
actions would be the mean field approximation, in which 
dimcnsionlcss magnetization per electron, fh = m/n/xs, 
is given by 



fh = tanh(6 + fh/t) 



(2) 



where t = T/T c , and density-dependent parameter T c > 
characterizes the interaction strength. We shall con- 
sider only T > T c , since this model implies a ferromag- 
netic instability at T = T c , forbidden in a 2D system 
with continuous symmetry. 

We solved Eq.[2] numerically to get fh as a function of 
b and t. The value of dM/dn can be expressed through 
fh and depends on an additional parameter, dT c /dn: 



1 dM 
Hb dn 



dm _ ~dmdT c 

n ^z = m - nt ^r^r ( 3 ) 



dn 



dt dn 



In the insulating phase (low densities) the interac- 
tion strength T c is expected to rise with n from zero at 
n = 0, due to increasing overlap between the localized 
electrons' wave functions. In the metallic phase (high 
densities) T c should decrease with density due to screen- 
ing. Positive dT c /dn is precisely the condition required 
to reproduce the peak in dM / dn(b). As an example 
we choose T c oc n k ; the results are qualitatively simi- 
lar for k = 1 and 2, see Fig. [4] We stress that the mean 
field dM/dn qualitatively captures all the features seen 
in Fig. la, including seemingly strange dM/dn > ug, 
and thus supports the claim of ferromagnetic interaction 
between spins in our system. 
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